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In this paper, we discuss the internal and external metric of the semi-realistic stars in relativistic 
MOND theory. We show the Oppenheimer-Volkoff equation in relativistic MOND theory and get 
the metric and pressure inside the stars to order of post-Newtonian corrections. We study the 
features of motion around the static, spherically symmetric stars by Hamilton- Jacobi mothod, and 
find there are only some small corrections in relativistic MOND theory. 
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• 1. Introduction 

, The "missing mass" problem comes of the discrepancy between two methods which relate to the measurements of 
mass to luminosity ratio. One is to measure the cluster's mass to luminosity by dynamic method; and the other is to 
directly measure the mass to luminosity through the rotation curve of spiral galaxies. To solve the problem, a prevalent 
way is to assume that there is much more invisible non-baryonic dark matter than visible baryonic matter in the 
universe. However, as an alternative to dark matter, Milgrom 1] proposed a Modified Newtonian Dynamics(MOND) 
which implies that Poisson equation should be rewritten as V (V$| V$|/a ) = 4nGp when |V$| <C a , where a « 
1 x 10 -8 cms~ 2 from empirical data. One found that MOND can explain not only the rotation curveQ of spiral 
£S| i galaxies but also the Tully-Fisher lawQ: Lk oc v*, where Lk is the infrared luminosity of a disk galaxy and v a is the 
^ ' asymptotic rotational velocity. 

Though MOND can explain the experiment successfully, it is not a theory and there exists some theoretical problems, 
such as the conservation of momentum and angular momentum. Recently, BekensteinQ has developed MOND to be 
\ a relativistic theory of gravity, abbreviated TeVeS, which depends on a tensor field g a p, a vector field ii a , a dynamical 
scalar field cf> and non-dynamical scalar field a. The theory also involves a free function F, a length scale I and two 
positive dimensionless parameters, k and K. Bekenstein has shown that TeVeS has general relativity(GR) as its limit 
when k — » with I oc fc~ 3 / 2 and K oc k in Friedmann-Robertson-Walker(FRW) cosmology, and TeVeS has a MOND 
and Newtonian limit under the proper circumstance. TeVeS also passes the elementary solar system tests of gravity 
theory |j| . Even though TeVeS is not a perfect theory, such as the violation of the local Lorentz invariance and the 
free choosing of the function F, TeVeS is still strong interesting recently. A recent paper discussed the large scale 
struture in TeVeS theory which show that it may be possible to reproduce observations of the cosmic microwave 
• • ■ background and galaxy distributions. 
. £h ' The external spacetime of the spherically symmetric objects in TeVeS has been explored by BekensteinQ and 
Giannios[p. Bekenstein has studied the post-Newtonian corrections of the static, spherically symmetric metric outside 
the spherically symmetric objects, and Giannios has derived the analytic expression for the physical metric. Up to 
now, there is no research on the stellar structure in TeVeS. In this paper we study analytically a semi-realistic star 
to get the primary characteristic of stellar structure in TeVeS. Firstly, we discuss the asymptotic vacuum spacetimes 
in the exterior of the stars by use of a static, spherically symmetric metric which can be transformed into the metric 
Bekenstein adopted. Secondly, we consider the Oppenheimer-Volkoff equation in TeVeS, and use the series expansion 
to compute the metric and pressure inside the stars. Finally, to get the observable correct effect of TeVeS, we explore 
the features of motion around the spherically symmetric stars by Hamilton- Jacobi mothod. 

2. The Brief Review of TeVeS 

The basic equations of TeVeS can be derived from the action 

S = S g + S v + S s + S m , (1) 
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where the actions of the tensor field(6 9 ), the vector field(6„), the scalar field(6,,) and the matter(S! m ) are: 



(2) 



6.1, 



6 S = ~ / [a 2 h a ^, a (/>, + -Gr 2 o^F(kGo 2 )\ (-gf^d^x, 



K 



32ttG 



(3) 
(4) 
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£(g^J a J Q ^,---)(-g) 1/2 d 4 x. 
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In the equations above, g is the determinant of Einstein metric g a p, R a p is the Ricci tensor of g a p just as in GR, 
h a P — g a/3 — iFif 3 , and g a p is the physical metric which is related to the Einstein metric by g a p = e^ 2 ^(g a p+il a ilp) — 
e 2 ^it Q ilg. F{p) is a free dimensionless function, with /i = kGcr 2 . A is a spacetime dependent Lagrange multiplier and 
f a symbolically represents the matter fields. The symbol | denotes the covariant derivatives with respect to g^ v and 
a pair of indices surrounded by brackets stands for antisymmetrization, i.e. A^B U ^ = A^B U — A V B^. 
Varying the action 6 with respect to g a ® , il Q and <$>, one will get the metric equation 



G. 



a/3 
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f afi + (i-e-**)Uf l f l 4 a iifi)+T a p +e Q/3 , 



the vector equation 



iOl [Q;/3I ;/3 + AiT + ^Ga 2 ^(j) t pg a ^(j) tl = 8ttG(1 - e^g^i&ftf, 



and the scalar equation 



[^<f> <a ], = kG[g aP + (1 + e-^iTil^ 



where 



T a /3 = O' 



Gr 2 a 4 F(fi)g a p, 



(6) 
(7) 
(8) 
(9) 
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and a pair of indices surrounded by parenthesis stands for symmetrization, i.e. A^B V ) = A^B U + A V B^. If one 
models the matter as a perfect fluid, then the energy-momentum tensor has the form 



T a p = pu a U/s + p(g a p 



yUp) 



(11) 



where p is the proper energy density, p the pressure and u a the 4- velocity, with u a — il Q , all three expressed in the 
physical metric. If varying the action S s with respect to a, one can arrive at the equation, —p,F(p) — i p 2 F'(p) = 
k£ 2 h al3 <f)^ a (j> y p. There is some discussion about the free function F(fi) in Because in this paper, we are interested 

in the post-Newtonian corrections not the purely MOND case, we shall take fi=l, i.e. a 2 — l/(kG). 

3. External Spacetime of a Static, Spherically Symmetric Star in TeVeS 

The metric of the static and spherically symmetric stars can be written as 

' )a pdx a dx si = -e v{r) dt 2 + e x{r) dr 2 + r 2 (d6 2 + sin 2 9dcp 2 ), 



ds A 



(12) 



where v and A are the function of r only. It is worth noting that metric Q12f) has the different form with that in 
Ref.0, but both metric can be transformed between each other. As we are looking for the static solutions, we will 
take the vector field to be pointing in the timelike direction. Normalized to il Q il a = —1, it has 



U a = {e-" /2 , 0,0,0}. 



(13) 
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Using Eqs.© and l|13|l . we get easily the components of gravitational equation. The tt component is: 

- 1 + + rAO = + ^^(^') 2 ] + * e"* + y - ^ + j] , (14) 



the rr component is: 



£(1 - e* + n/) = 8nG[pe^ + -^\V) 2 ] + Ke~\- (15) 



and the 99 component is: 

The ipip component is proportional to the 99 component, so there is no need to consider it separately. In principle, 
the first two components of gravitational equations are enough to solve the metric. Furthermore the scalar equation 
becomes 



e ^]' = kG(p + me^e-^. (17) 



Similar to Bekenstein's definition^, we define the "scalar mass" m s as 

r R 

m s = m s (R) = 4tt / (p + 3p)e (,y+A)/2 e- 2 Vdr, (18) 
Jo 

where R is the radius of a spherical symmetric star. 

We consider the case that the exterior of the star is vacuum, and have p = and p — 0. For r > R, Eq. l|17|l becomes 

0' = fc , s6 2 ■ (19) 



Substitute p, p and <f>' into Eas. (|14|) and 1)150. we have 



a lft ,(Gm,fe-" r ^ Ar (i/) 2 f'A' *A 



and 

— (1 - + „/) = fc^^ — + if e- A [ - . ( 2] ! 



For getting the post-Newtonian corrections, one can expand e v and e as 



e 



" = a [l--+a 2 (-) +«3(-r + -.], (22) 



e A = /3 [1 + /?!(^) + /? 2 (^) 2 + /? 3 (-) 3 + ■ • ■] , (23) 



where r e is a length scale to be determined in the next section, and the size of the coefficient of the r e /r term in 
Ea. (|22|l has been absorbed into r e . Substituting them into Eas. (|20|l and (|21() . matching the coefficient of like powers 
of 1/r, and solving the recurrence relation, we obtain 

a 2 = 0, (24) 
if k (Gm s ) 2 1 



Q!3 



48 a Q 4tt 6r 2 ' " 
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= 1, 
= 1, 

= 1 



(25) 



A' fc (Gto s ) 2 1 

ir ~ ~ 4tt ^f' 

5 A fc (Gm s ) 2 5 
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To gain the value of the coefficient of ao, we can use the fact that the limit of TeVeS is GR when k — > and K — > 0. 
Comparing the expansion with the Schwarzschild metric in GR, we get — 1. After substituting the expansion into 
Ea.p7[l and integrating the equation, we have 

kGm s 1 kGm s r e 1 _ 3 

<P = <Pc -j 3 2+°( r )' ( 26 ) 

47T r 57T r z 

where the constant <j) c is a cosmological value and can be absorbed in the rescaling t and r coordinates. 

Because of the relation g u = —e 2 ^ +u , g rr — e _2 ^+ A and ggg = g vv / sin 2 9 — r 2 e~ 2 ^, the physical metric becomes 

9u = -1 + ^ - 2G WGW™~fr) + (r-% (27) 

~ 9rr = i | 2G ^ m | ZGW^m+^-fc^ + frg | 

rr r r 2 ' 

5ee = r 2 + (2G A rm-r e )r + 2G J vm(G A rm-^)+0(r- 1 ), 

g w = sin 2 Ogge, 

where 2Gjvm = r e + 2fc^. 

4. Oppenheimer-Volkoff equation in TeVeS 

To get some characteristic of stellar structure in TeVeS, we will research a semi- realistic stars in this section. In 
other words we can get the kernel information about the stellar structure in TcVcS through the simplified model, and 
to find how difference the stellar structure is between in TeVeS and in GR. The simplified model of star comes from 
assuming that the fluid is incompressible: the density is a constant p~o at r < R and it vanishes at r > R. In the 
internal of a star with ideal fluid T a f3, the r component of energy- momentum conservation law T a ^ — gives 

2 Tr- (28) 

Rewriting Eq. ltT4"|) and Eq.JTSJl, we have 

-A 



— ( - 1 + e A + rX') = 8nGj5, (29) 



and 

-(l-e x + rv') =8irGp. (30) 



where p = g#* + sf^W + £ - + and p = [pe~ 2 * + ^{<t>>?] + ^e">[ - . 

According to Ea. (|29|l . it is convient to replace A(r) with a new function m(r), given by 



or equivalently 



Then Ea. l(2§J) becomes 



which can be integrated to obtain 



m(r) = i (r - re A ) , (31) 



e-^!- 2 ^). (32) 



dm(r) _ , 



dr 



4irr 2 p, (33) 



m(r) = 4vr / p(r')r u dr'. (34) 
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For r>R, 



e- x = 1 



2G 



An f* p(r')r' 2 dr' + 4tt f R p(r')r' 2 dr' 



1 _ 2Gm g (R) _ 87rGf R p(r'y 2 dr' 



where m g is the "gravitational mass", with ra g = m g (R) = Ait J Q p(r')r' 2 dr' , and in the integration the expression of 
p{r') depends on the internal metric which has a complex form. At r — R, using Eas. (|15|) . ((25(1 and 1(35|) . we have the 
relation between m s , m g and r e 

Kr 2 kG 2 m 2 /r e 3 \ 

re - ^iT + -SS- + ° (^J = 2Gm - (36) 

In terms of m(r), Ea. (|30l) can be rewritten 

dv _ 2Gm{r) + 8irGpr 3 
dr r [r — 2Gm(r)] 

Combining Eq. ((37(1 with Eq. ((28(1 allows us to obtain the Oppenheimer-Volkoff equation 9] in TeVeS 



(37) 



dr ~ [pe +pe > r[r-2Gm(r)} ' [M) 
To post- Newtonian corrections, in the interiors of static, spherical symmetry stars, we expand e", e A , <j> and p as 

^ = «o [1 + (f ) 2 + « 2 (f ) 4 + « 3 (f ) 6 + - - ■] , (39) 

e A = 6q [! + 6l (L) 2 + 6 2 (^1) 4 + b 3 (Lf + ...], (40) 

7*2 

^ C o[l + Cl (f) 2 +c 2 (f) 4 + c 3 (f) 6 + ...], (41) 

^rfi(f) 2 + d 2 (f) 4 + d 3 (f )%..., (42) 

where r,- is a length scale to be determined later and the size of the coefficient of the r/rj term in Eq. 139(1 has been 
absorbed into r t . In Ea. (|42|l . we need not write down the constant term because the constant is a cosmological value 
and can be absorbed in the rescaling t and r coordinates. Substituting them into Ea. l|14[) . Ea. 1(15(1 . Ea. 1(17(1 and 
Eq. ((28(1 , matching the coefficient of like powers of r 2 , and solving the recurrence relation, we have 

_27k(2-K) 8ttC?po(5 - AK)rj 5 - AK 2 + 6kGp r 2 + K(8 - 3kGp r 2 ) 
a ' 2 ~ 80^ + 30(2 -K) + 10(2 - K) ' { ' 



b, = 



1, (44) 

SnGpo 2 
h = K + — - — 1\ , 

K 3 (16tt - 9k) K 2 (-3A8n + 243fc + 6A0n 2 Gr 2 p o ) 
16ir{~2 + K) + 120tt(-2 + K) 

K[-2A3k - 16327T 2 Gr, 2 / 5 + I280ir 3 G 2 rfp 2 l + 72^(2 + 3fcGr 2 p )] 12807T 3 G 2 pg + 27fc(3 + 8irGr 2 p ) 
180tt(-2 + K) 90tt(-2 + K) ' 

3 snGrf 
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1 kG(2-K) (10 - 5K)[8nG - 3kG(2 - K)] 
3 + 16ttG + 48ttG(-2 + K + 8nGp rf) ' 

27k 2 (-2 + K) a 32{-5 + 4K)ir 2 G 2 rfp 2 

128tt 2 (-6 + 3K + HirGr^po) + 15(-2 + K)(-6 + 3K + 8irGrfp ) 
47rGrfpo[-15 + A8kwGrfp + K(27 - UkGrfpp) + 2K 2 (-3 + 5kGrfp )] 

15(-2 + K)(-6 + '3K + 8TrGr 2 p ) 
3k(-2 + K) 2 (5K - 2(9 + bkGr 2 p )) 

40tt(-6 + 3K + 8TrGr 2 p ) 
60 - 492kGr 2 p Q + 40k 2 G 2 rfp 2 + 8K 2 (3 + WkGrfpp) + if (-93 + 86kGr 2 p - 20k 2 G 2 rfp 2 ) 

40(-6 + 3A + 87rGr?p ) 



di 

rf 2 



3Ar(2 — if) 



8tt 



(46) 



1927T 



[-12isT 2 + 6(-3 + 2(4tt - k)Grfp ) + K(33 + 2(-16tt + 3k)Grfpo)}, 



To get the value of the coefficient of ao, we also use the fact that the limit of TeVeS is GR when k — ► and 
K — > 0. Comparing the expansion with the internal metric of a static, spherically symmetry star in GR we have 

5 9Gm, 



ao 



2R 



1 



K 



irGp a r 2 
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At r = R, using Eq.JTSJ, Ea. (|4^|) and Ea. l|33|l . we have the relation between m g and 



(47) 



18K 2 (4ir - 9k) 4K[-27k + 32ir 2 Gp r 2 + 8tt(2 + 3kGp r 2 )] 24fc(3 + 8TrGp r?) \ R 



n{-2 + K) 80tt(-2 + A) 
2Gm a . 



80tt(-2 + K) 



80tt(-2 + if) 



Using the relation g tt 



(Jit 



-e 2,f > +u , g rr = e - 2 '> }+x and g ee 



9 W I sh 



re 



-20 



, we have 



9Gr 



2R 

9Gm g 
2R ' 



3 I 2Gm Q 
2 Jl - 



R 



(48) 



1 - 



2Gm 



K 



3k(-2 + K) 8 
8^ + 3 

3k(-2 + K) 



R 

nGpoT 2 



1 - 



3fc(-2 + A) 



8- 



O 



O 



8tt 



-] +0 

Ti 



9<p<p 



-.2 fl- 



it is clearly that the internal physical metric of spherical star in TeVeS will get back that in GR when k — ► and 

Though the stellar model discussed in this section is not a perfect model, it is still a first-order approximation of 
small star of which the pressure is not large. The procedure above shows us that theoretically we can deal with a real 
star in TeVeS in the similar way, but the numerical calculation will be more complex. 



5. Features of Motion Around a Spherically Symmetric Object in TeVeS 



Based on the analyses of the sections above, the features of motion around any spherically symmetric compacted 
objects not only the black holes can be discussed now. An analytic solutionjjj has been gotten in spherically symmetric 



7 



spacetime in TeVeS for arbitrary k and K, where the metric in TeVeS is consistent with these in general relativity, if 
k = and K = 0. The solution is based on the line element 

ds 2 = g a/3 dx a dx f} = -e u{p) dt 2 + e^[dr 2 + p 2 (d9 2 + sin 2 9dip 2 )}, (49) 

which is equivalent to the line element in Ea. (|12fl . Between Ea. (|12fl and Ea. (|49|l . it is easy to find the relation: 

e »( P ) = e «'(r) j ( 50 ) 

e^dp 2 = e A «dr 2 , 
e^p 2 = r 2 . 

And the static, spherically symmetric analytic solutions are: 

W^c + ^lnf^V (51) 



%itp c \p + Pc 

<°-{hr) • (52> 

P \P + PcJ 

where p c = ^f^Jl + ^(^p^) 2 — "f- Pa * s a l en gth scale to be determined. And m s is the "scalar" mass which is 

defined by m s = An J^(p + 3p)e l/ / 2+3? / 2 ~ 2< ^ p 2 dp, where R is radius of the matter's boundary^]. Using the results of 
post-Newtonian corrections, one has the relation between m s ,m g and p g 

^- 3 -ji- k -^r +0 ^ /R2) - 2Gm - (54) 

where m g is the "gravitational mass". For the sun p g /R ~ Gm s /R ~ 10 , one will find that p g w 2Gm g with 
fractional accuracy much better than 10 -5 . The physical metric may reduce to gu = gtte 2 ^ and gu = gue~ 2 ^. And 
one will arrive at 



~9u = ~ "—^ (55) 

\P + Pc' 



~9 PP = ^ = ^= [A ^ L ( E ^) (56) 
p z p 2 sir 9 p 4 \p + p c 1 



where a = + fc 4 ^™ a . When k = and K = 0, the metric will reduce to the metric in GR. 

The relation between energy and momentum for a test particle of rest mass m in curved space is 

~g a(3 p a pp + m 2 = 0, (57) 

where p a = dS/ dx a and S is the Hamilton- Jacobi function. Thus Hamilton- Jacobi equation for propagation of wave 
crests in TeVeS spherically symmetric geometry becomes 

3"(f ) 2 +^(g) 2 +9 PP p- 2 (§) 2 +9 PP P- 2 s™ 2 ^) 2 rn 2 = 0, (58) 

where p = ^/g rr r. Set 

S=-Et + S 1 (p) + S 2 (6)+S 3 (<p). (59) 

Substitute Ea. (|59"|) into Eq.JSBJ, we have 

~ g u& + cf(^) 2 + ~9 pp P- 2 {^) 2 + 9 PP P- 2 sin" 2 9(^) 2 + m 2 = 0. (60) 
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Solving the equation we find the solution of Hamilton-Jacobi function 



S = -Et 




fJ 



pp 



-E 2 



— dp 



L 2 



sin 2 ( 



■d6 



P<p<P, 



(61) 



where L is the angular momentum and E is the energy of the system |8j . 

Find the relation between t and p by considering " interference of wave crests" belonging to slight different E values 



Then one has 



dS 

—^ = -t + 
dE 



dp 
~dt 



2—E 

g pp 



_9^_p2 _ rn 2 . _ L 2 

gPP gPP ~~p2 



-.dp = 0. 



(62) 



.911 E 2 

g pp 



L 2 



-M^E 
g pp 



(63) 



When dp/dt = 0, Eq. I|63|) infers the effective potential energy of the system: 

1/2 

cipp I mf Ij" \ 

U = 



gPP 


' m 2 




gtt 







(64) 



Substituting the physical metric g a and g pp into Ea. H64() . we have the effective potential for unit mass 



U = 



Pc 



P + Pc 



a/2 



1 



P 



P- Pc 



171 (P 2 -Pc) \P + Pc 



(65) 



where U — U jm is the effective potential for unit mass of the test particle. Defining three dimensionless quantities 
t = p/p g , v c = p c /pg and L = L 2 / (m 2 p 2 g ), we get 



where a = 



2t c ( 



1 



k Gm £ 



u 



and r c = 4 



a/2 



L 2 - 



1/2 



(66) 



1 ' 7T V Pa ) 2 



Energy, in units of the rest mass m of the particle, is denoted E — E/m. It is worth noticing that one can define a 
turning point by the condition U 2 = E 2 . Obviously, the minimum and maximum of the potential depend on the L. 
The roots of dll/dp = are given in terms of the L by numerical calculation, so that we can define the lowest reduced 
angular momentum L cr a by the condition U m i n (Li owest ) — U max (Li owest ). In other words, there is no periastron for 
L < Liowest, and any incoming particle is necessarily pulled into p — p c . The radius of stable circular orbits represents 
the point sitting at minimum of effective potential. The last stable circular orbit corresponds to Li owest . It is easily 
found that Li owest and the radius of stable circular orbits depend on k and K . We also define L cr a by the condition 
U m ax{L C rit) = 1, then there are bound orbits for L > L cr u and test particles coming in from p — oo with E 2 < U 2 lax 
reach periastrons and then return to p = oo, but particles from r = oo with E 2 > U 2 iax get pulled into p = p c . There 
are unstable circular orbits at the maximum of the effective potential, in which maximum moves outward form the 
radius of stable circular orbits for L — oo to the radius of stable circular orbits for Li owest . 

Giannios 7] had proved that the solution for the metric of static, spherically symmetric black holes is identical to 
the Schwarzshild solution in GR. To show the difference of features of motion around a spherically symmetric object 
between in TeVeS and in GR, we can discuss the corrections of stars, such as sun, but not black holes in TeVeS. For 
sun Gm s /pg ~ 1O~ 1O 0, so it is reasonable to neglect the terms with Gm s /p g in the expression of a and r c , i.e. k is 
insensitive to Li owest and L cr i t . We define the relative errors of Li owest and L cr i t 



AL 



lowest 



AL r 



Llowest{K) — Li owest (K — 0) 
Llowest(K = 0) 

= L crlt {K) - L crlt {K = 0) 
L crlt (K = 0) 



(67) 



(68) 
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TABLE I: Values of Li oweat , ALi owest and L cr u, AL cr i t for different K 



K 


■^lowest 






Lcrit 







l 


.73205 







2 





0.01 


l 


.73221 





.009013% 


2.0002 0.00992684% 


0.02 


l 


.73236 





.018022% 


2.0004 


0.019847% 


0.03 


l 


.73252 





.027025% 


2.0006 


0.0297605% 


0.04 


l 


.73267 





.036023% 


2.00079 


0.0396673% 


0.05 


l 


.73283 





.045016% 


2.00099 


0.0495674% 


0.06 


l 


.73299 





.054003% 


2.00119 


0.059461% 


0.07 


l 


.73314 





.062985% 


2.00139 


0.0693478% 


0.08 


l 


.73330 





.071961% 


2.00158 


0.0792281% 


0.09 


l 


.73345 





.080932% 


2.00178 


0.0891018% 


0.1 


l 


.73361 





.089899% 


2.00198 


0.0989688% 


0.2 


l 


.73516 





.179271% 


2.00395 


0.19728% 


0.3 


l 


.73669 





.268126% 


2.0059 


0.292944% 


0.4 


l 


.73822 





.356470% 


2.00784 


0.391974% 


0.5 


l 


.73975 





.444313% 


2.00977 


0.488381% 


0.6 


l 


.74126 





.531661% 


2.01168 


0.584175% 


0.7 


l 


.74276 





.618522% 


2.01359 


0.679367% 


0.8 


l 


.74426 





.704904% 


2.01548 


0.773968% 


0.9 


l 


.74575 





.790813% 


2.01736 


0.867986% 


1.0 


l 


.74723 





.876256% 


2.01923 


0.961433% 


1.1 


l 


.74870 





.961241% 


2.02109 


1.05432% 


1.2 


l 


.75016 


1 


.045774% 


2.02293 


1.14665% 


1.3 


l 


.75162 


1 


.129861% 


2.02477 


1.23843% 


1.4 


l 


.75307 


1 


.213508% 


2.02659 


1.32968% 


1.5 


l 


.75451 


1 


.296722% 


2.02841 


1.4204% 


1.6 


l 


.75594 


1 


.379509% 


2.03021 


1.5106% 


1.7 


l 


.75737 


1 


.461874% 


2.03201 


1.6003% 


1.8 


l 


.75879 


1 


.543824% 


2.03379 


1.68948% 


1.9 


l 


.76020 


1 


.625363% 


2.03556 


1.77817% 



With the numerical calculation, the values of Li owes t, ALi owes t and L cr it, AL cr u for different K are listed in TablcQ] 
Note that K = corresponds to GR. We find Li owest and L crit have a small difference between GR and TeVeS. 

6. Conclusion 

Based on relativistic MOND theory, we have studied the spacetime in the external and internal of a semi-realistic 
star with a constant density inside the star in this paper. We show the Oppenheimer-Volkoff equation in TeVeS, and 
calculate the metric coefficient and pressure inside the star to post-Newtonian order. Therefore, one can determine 
the relation between m s , m g and R by integrating the metric and scalar equation inside the star. The more complex 
equation of state to approach the stellar structure is left for future work. Furthermore, we study the features of a 
test particle around the star in TeVeS and get the effective potential which can determine the features of particle's 
motion. We define two relative errors to discuss the corrections of TeVeS and show that there is a small correction of 
TeVeS which can be tested by the future experiments. 
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